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We present a unified account of magnetic exchange and 
Raman scattering in the quasi-one-dimensional transition- 
metal oxide NaV2 05. Based on a cluster-model approach 
explicit expressions for the exchange integral and the Raman- 
operator are given. It is demonstrated that a combination of 
the electronic-structure and the Dzyaloshinskii-Moriya inter- 
action, allowed by symmetry in this material, are responsible 
for the finite Raman cross-section giving rise to both, onc- 
and two-magnon scattering amplitudes. 

PACS numbers: 75.30.Gw, 75.10.Jm, 78.30.-j 



I. INTRODUCTION 

More- than four decades ago, DzyaloshinskiiS and 
Moriyaa showed that the inclusion of spin-orbit coupling 
into the description of low-symmetry magnetic systems 
generates an anisotropic exchange interaction, the so- 
called Dzyaloshinskii-Moriya (DM) interaction. 

In the early nineties, this interaction was discussed 
intensively in connection with the copper-oxide com- 
pounds. In particular, La2Cu04 exhibits a small gap 
in the spin-wave spectrum and a finite net ferromag- 
netic moment in each plane due to an out-of-plane cant- 
ing of the spina. These features were attributed to DM 
interactions.H13 Yildirim et al.a did a careful microscopic 
study of this mechanism for tetragonal copper-oxide sys- 
tems. In particular, their analysis proved that the or- 
thorhombic distortion present in these materials is irrele- 
vant to the out-of-plane magnetic anisotropy. Moreover, 
they showed that not only the antisymmetric anisotropic 
superexchange between two neighboring spins is impor- 
tant but the symmetric one as wellHia 

The DM interaction has gained renewed interest in 
the context of the novel transition-metal oxide NaV205, 
which is believed to be a quarter-filled ladder compound 
in its high-temperature phase .0 At Tq = 34K a phase 
transition, the interpretation of which is still controver- 
sial, takes place in this material where charge ordering 
(2V +4 ' S — * V +A + V +5 ) occurs simultaneously with the 
opening of a spin-gap of approximately 10 meVJij A se- 
ries of recent studies has addressed the nature of the low- 
temperature stateE3c3. 

In this context it is of interest that_s£exjy recent Electron 
Spin Resonance (ESR) experimentsliatfl have detected a 
considerable anisotropy of the absorption intensity with 
respect to the magnetic field orientation, which has been 



attributed to the DM interaction. Apart from ESR, Ra- 
man scattering in the presence of a magnetic field is an 
alternative experiment for the observation of possible ef- 
fects due to DM interactions. Unfortunately however, at 
present, the various experimental settings in search for 
such effects ia, Raman scattering have been unsuccessful 
for NaV2 05.EJ This may not be conclusive yet, since the 
Raman cross section of the relevant scattering process 
could be too small to be observable. 

In order to shed some light onto this scene, a mi- 
croscopic analysis of the magnetic exchange and Ra- 
man scattering-operator seems highly desirable. How- 
ever, apart from-early work specific to the copper-oxides 
superconductorsEil such analysis is lacking. Therefore, it 
is the purpose of this paper to present a detailed descrip- 
tion of the Raman operator for NaA^Ov In this context 
particular emphasis will be given to the role of the DM 
interaction which by symmetry is allowed in this mate- 
rial. 




FIG. 1. Crystal structure of NaV2 0s in the high-tempera- 
ture phase. The star denotes the location of the center of 
inversion, the dashed lines the constituting V-O-V ladders. 
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II. HAMILTONIAN 

Above the spin-charge transition temperature Tc 
NaV20,s-<;rystallizes in the centrosymmetric Pmmn space 
group.BtJ The compound consists of VO5 square pyra- 
mids sharing edges in the ab layer and chains of Na lo- 
cated between the ab layers. The superexchange inter- 
action between vanadium sites is mediated through the 
pyramid's base oxygens and the relevant structural ele- 
ment of NaV2C>5 can be thought of consisting of ladders 
of V-O-V rungs along b which are weakly coupled along 
a (see Fig. |l|). 

Discarding single-ion anisotropy, a general form of any 
scalar two-spin interaction between consecutive rungs 
along the 6-direction of the ladder in NaV20s consists 
of two contributions 



= JH {Heis) +DH {DM) 



(1) 



i.e. the isotropic Heisenberg exchange H^ Hels ^ and the 
Dzyaloshinskii-Moriya interaction H^ DM >: 



I 

ff(™) = ^e.(S,xS l+1 ) 



(2) 



where Si denotes the total spin on rung /. The form of the 
Dzyaloshinskii-Moriya (DM) vector e is determined by 
requiring that the energy of any configuration of spins has 
to be invariant under the symmetry transformations of 
the crystal structure. In our case, crystallography allows 
for a DM vector along c, i.e. e = (0,0,1). Note that 
this vector is defined locally in each unit cell and is not 
forbidden by the inversion center of the crystal structure, 
which lies in between two V-O-V ladders (see Fig. |l|). 

The standard derivation of (||) for the case of magnetic 
moments localized at single ionic sites can be found in the 
literatureB'Lrlj where it is shown that both terms in this 
Hamiltonian can be derived from generalized exchange 
processes. In the present paper, however, we aim at a 
microscopic derivation of (^) for the case of the mixed- 
valence system NaV20s. We start from a three-band 
Hubbard- model for NaV2C>5 in which we retain only the 
two active ci-orbitals of the V- and the single 0(l)-site on 
each rung. For simplicity we consider the 0(2)-sites on 
the legs of the ladder to be integrated out, giving rise to 
an effective V-V hopping —t\\ along the leg. We denote by 
d, the creation operators for spin-a electrons in V-d xy 

orbitals on site a (a = 1, 2) of the Vth rung and by p\ 
the creation operator of spin-er electrons in the 0(\)-p v 
orbital on the Z'th rung. The Hamiltonian H = Ho + Hi 
reads: 



La 



Hi = ( t„ d) d,,-, + t* d], , d, 

1 / j \ ° I. a cr l-\-la<T cr l-\-laa lac 
l,acr 

X] (PlvPl+lv + Pl+laPla) ' ( 3 ) 



<-pp 



The hybridization matrix elements in Eq. ([}]) are: (i) 
t± which denotes the direct hopping of electrons between 
the ~V-d xy orbitals on sites 1 and 2 on a rung, (ii) td P 
which denotes the charge transfer integral between the V- 
d xy and 0(l)-p y orbitals on a rung, (iii) t a = — 1\\ + i<j\ is 
the hopping of electrons with spin a — ±1 between the V- 



d xy orbitals on two consecutive rungs along the ladder di- 
rection b. The spin-dependence of this hopping-integral is 
allowed due to the lack of a center of inversion in-between 
two rungs and arises from the spin-orbit couplingcl of 
strength A, where A ~ A. The transfer matrix-elements 
t a are diagonal in the spin quantum numbers because we 
have chosen the quantization axis for the spin to be along 
c, i.e. the main crystallographic axis, and finally (iv) t pp 
denotes the hopping of electrons between 0(l)-p y orbitals 
on two consecutive rungs along b. The spin-dependence 
of tpp is small and will be discarded in the remainder of 
this paper. _ 
The parameters involved in (||) have been estimatedtl 



to be e % 



-3eV, t 



0.25 eV, t dp 



-leV, U 



-0.175 eV, tpp w 0.5 eV and U « 2.8 eV. The Coulomb 
repulsion U leads to the formation of local moments on 
the rungs and in the following we will study the inter- 
action between these local moments. For simplicity we 
consider the case U — > 00, since we expect that any finite 
U will lead to qualitatively similar results whila4ncreas- 
ing the complexity of the calculation needlesslylilJ 



III. EXCHANGE COUPLINGS 

As a first step towards the evaluation of the exchange 
couplings J and D, we diagonalize the hamiltonian Hq 
for an isolated rung. The ground state belongs to the 3- 
particle subspace. In addition, for the calculation of the 
exchange matrix elements, intermediate states in the 2- 
and 4-particle sector on a rung are required. Details of 
the derivation of the relevant eigenstates and eigenvalues 
of this cluster problem are stated explicitly in Appendix 
A. The exchange matrix elements are obtained by con- 
sidering the process which describes a spin-flip between 
two consecutive rungs up to second order perturbation 
theory in H±, as defined by J~ + and J H in Eq. (|AS|). 
In particular, 



+ ,/+-) = Re(J-+). 



(4) 



Substituting the electronic m odel- param eters cited in the 
previous section into Eqns. (A13) and (A15) we obtain 
J w 0.049ey w 568J£] which agrees very well with the 
experimental valueOO of « 560A". 
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Noting that t± — tZ and that D ■ (Sj x S;+i) 
(D/2i) (SfSf +1 - SfSf +1 ), we have 



Im{J~ 



(5) 



for the DM- coupling. MoriyaB has estimated that the 
order of magnitude of D should be 



D ~ (Ag/g)J 



(6) 



where g is the gyromagnetic ratio of the Vanadium ion in 
octahedral crystal symmetry and Ag is the correspond- 
ing deviation from the free-electron value. By considering 
the g-values obtained from ESR measurements!!^, we ar- 
rive at Ag/g sa 0.01. Then, from Eq. (g), Eq. (§) and 



Eq. ( A15) where we evaluate Dasa function of A, we get 
an estimate for A, i.e. A « lmeV. 



IV. RAMAN SCATTERING 

Fleury and LoudorH have shown that light scatter- 
ing from a spin system, depending on the polarization 
geometry of the incoming and outgoing electric fields, 
can lead to inelastic photon-induced supercxchange. This 
has established Raman scattering as an important probe 
to obtain information on the local exchange dynamics 
in magnetic systems complementary to inelastic neutron 
scattering (INS). In the following we will generalize the 
early ideas of Fleury and Loudon to the case of NaV2 0s 
clarifying the role of the DM interaction. In particular we 
find that in the case of a polarization of both, the incom- 
ing and outgoing photon fields parallel to b, i.e. along the 
legs of the ladder, the Raman scattering operator H^- R ' 
can be expressed as 



+D R (u; m , out )H^ DM ^ 



(7) 



where the Ui n and io ut are the frequencies of the incom- 
ing and outgoing photons. The microscopic derivation 
H^ R \uJin iOU t) is placed into the appendix. It is identical 
to that of the magnetic exchange integral with however 
the virtual hopping into the intermediate state of the 
exchange process driven by the coupling of thc-szector 
potential A = (0, Af,, 0) to the current operatoio, i.e. 
H\ of (g) has to be replaced by j • A with the current 
operator j = (j a ,jb,jc), 



3b = *c 



PlaPl+la Pl+laPla 



E( 



The total magnetic Raman scattering ampl itude is then 
given up to second order in %je by Eq. ( A10 ) of Ap- 
pendix A. From this, the definition of J R (uJi n ,out) and 
D R {u) in , out ) is analogous to (f|) and (|) 



Jn(Uin,out) = Re(R + (LUir h out)) 

D R (u> in<(m t) = Im(RT + (uJin,out)) 



(8) 



Note that a magnetic Raman process is possible only if 
H^ R \u>i n ou t) induces transitions between different eigen- 
states of H ( S )@"R i.e. if 



[H( R \w ini0ut ),H^]-. 

-D r (uj 1 



out )J)[H^ s \H^]^Q (9) 



^From (^) we conclude that magnetic Raman scattering 
from NaV 2 5 , if modeled by (g), arises because two con- 
ditions are simultaneously satisfied. First, the existence 
of a spin-orbit coupling leads to a non- vanishing commu- 
tator in (||). Second, because the number of available 
paths for the magnetic and the photon induced exchange 
is larger than one and because Ui nou t ^ the factor 
of J R {LOin,out)D - D R (uj m , out )J is nonzero. The latter 
is true despite the formal similarity between the Raman 
scattering amplitude and the magnetic exchange integral, 
because H^ R \uji n ^ ou t) displays an additional dependence 
on the photon energies. More specifically, for a single 
exchange path J R (u>i n ,out)D - D R {uj in<out )J = for any 
value of LUin^out while for more than one exchange path 
jR{uin,out)D-D R (uj in , out )J vanishes only at uj m ^ out = 0. 

Next, we would like to point out that from the two 
terms, H^ Heis ^> and H (DM 1 which make up H^ R \u) int0ut ) 
it is actually ff( DM ) which drives the magnetic Raman 
process. Up to now we have only considered anisotropic 
contributions to fi^ to leading order in A. Kaplantl and 
Shckhtman et alB have shown that, in general, the next- 
order term H^ KSAE ^ ~ J2i( e ' s 0( e ' s *+i) contributes to 
the spin Hamiltonian with a very specific prefactor: 



H (s} = J H {Heis) +DH {DM) 
V J 2 -D 2 -j) H^ KSAE \ 



(10) 



Using this it is possible to transform H 1 -^ into an equiv- 
alent Hamiltonian of the plain Heisenberg form using the 
unitary mapping 



Sf= cos tpi Sf — sin (pi Sf 
Sf= - sin if 1 Sf + cos if 1 Sf 



(11) 



with Sf — Sf, tpi — 2lip , and tan(2<p ) = D/J. Ex- 
pressed in terms of S; the Hamiltonian reads H 1 --^ = 
\/ J 2 + D 2 J^i Sz • S;+i. Now, we note that higher 
order terms in A will also contribute to the Ra- 
man operator. However, following the discussion af- 
ter (g) it is obvious that ( |Tl| ) will not simultane- 
ously reduce and H^ R > (wj n ,o«t) to a canonical 
Heisenberg form. Therefore, in the new basis, the 
Raman operator takes on the form H^ R \u}i n}OU t) — 

j R (LO m ,out) H( He ^ +D R (UJ 

in, out 

)H(DM) +0 (\ 2 ). The 
only part of H^ R ' (cui nMUt ) which does not commute with 
to lowest order in A is the DM-interaction, i.e. 

H (fi> k,» t ) = D R (Uin,out)H (DM) - 
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This completes our derivation of the Raman opera- 
tor for the homogeneous phase of NaV^Os as realized 
for T > Tq. Quite generally the preceeding demon- 
strates that a DM-contribution to the Raman opera- 
tor ifw ((*>in,out) will occur in multiband systems when- 
ever a DM exchange-interaction is allowed locally. Ob- 
viously it is tempting to analyse the effects of this form 
of H^ R \u>in tOU t) also on a dimerized spin-liquid state, as 
present in NafcOs for T < Tc and similarly in CuGeC>3 
for T < T S pB. To this end let H {R \w in>ou t) act on 
a pure dimer state |$ ) = |si ...s M . ..) , where fi labels 
nearest-neighbor pairs of spins which are in a relative 
singlet state |s M ) - for the case of NaV^Os these pairs 
of spins correspond to pairs of rungs (21,21 + 1). One 
obtains 

ff(^)|$ )= £(-2i|...i*...)-|-^ + i-> 

+ ■ (12) 

Here \t") (a = x,y,z) refers to triplet states on the 
dimer-bonds. While the 2nd and 3rd term on the rhs. of 
( |l2| ) comprise of the usual total-spin zero, two-magnon 
excitation, the first term refers to a single-triplet state 
of only z-direction. This shows that single-magnon 
Raman-excitations are allowed in the presence of the 
DM-interaction. A single-magnon Raman line of this 
type has a clear experimental signature: it should show 
no splitting in an external magnetic field parallel to e 
(here along z) and it should split into two branches for a 
field perpendicular to the DM-vector. To our knowledge 
this signature has not yet been observed in experiment. 

V. CONCLUSIONS 

Motivated by recent ESR experiments00 which probe 
the existence of a DM interaction in NaV^Os, we have 
presented a microscopic study of the possible impact of 
this interaction on the Raman process. We have derived 
the Raman operator in the homogeneous state of NaV^Os 
and, additionally, have discussed its effect in the dimer- 
ized state. 

In the dimerized state two Raman-modes have been 
observed in NaV^Os in bb-polarization at 66cm -1 and 
104cm -1 . Tentatively these modes have been ascribed to 
magnetic bound states of total-spin zero.E2l On the other 
hand, for T < Tc INS displays two well defined magnon- 
excitations, the energies of which, if properly zone-folded 
to zero momentum coincide with the aforementioned two 
Raman modeal3'Ej. Yet, Raman experiments show no 
indication of a splitting of these modes in an external 
magnetic field. We therefore conclude that the Raman- 
modes should result from a two-magnon processes (see 

Clear evidence for a DM-vector in NaV^Os along the. 
z-direction has been provided by ESR experiments.t£l 
While these authors have interpreted there findings in 



terms of quasi-static charge-fluctuations above Tc, we 
believe, in view of the results presented here, that such 
an interpretation of the ESR-data is not necessary. In 
fact, the ESR-cxperimcnts can be understood in terms of 
the local DM-vector present also in the high-temperature 
phase. 

In conclusion we have pointed out, that a local DM- 
vector gives rise to a non-trivial DM-contribution to 
the magnetic Raman-process whenever at least two non- 
equivalent exchange paths exist between the two mag- 
netic moments considered. We have presented an ex- 
plicit evaluation of this DM-contribution to the Raman- 
operator for the case of the quarter-filled ladder com- 
pound NaV205 and we have shown, that one- and two- 
magnon processes arise naturally within this scenario. 
We have obtained estimates for the exchange-coupling 
constant along the 6-direction in good agreement with 
experiment. Moreover we have evaluated the spin-orbit 
coupling constant within our cluster-approach. Finally, 
we note that evidence for DM-interactions-in the two- 
dimensional dimer-compound SrCu2(B03)2L3rhave been 
found by ESRCj and far-infrared spectroscopyO. There- 
fore one might speculate if one-magnon Raman modes 
with the special signature described in the previous sec- 
tion could be observable in SrCu2(B03)2- 
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APPENDIX A: EXCHANGE INTEGRAL AND 
RAMAN AMPLITUDE 

In this appendix we present details of the evaluation of 
the 2, 3, and 4-particle eigenstates on a rung as well as 
the matrix elements relevant to the exchange integral and 
the Raman operator. We begin with the 3-particle space 
on rung I, which, in the subspace of no double-occupancy 
of the d- levels and total-spin z-component S z =| can be 
created by 

3' — rft rft 3' — r/l d^ 

= ^2llPf\Pli ^4/| = ^U^PltPll (Al) 

The set of corresponding states with S z =!, is obtained 
by reversing "f to J. for each operator without changing 
their relative order. Diagonalizing the rung-Hamiltonian 
in this sector yields a (spin degenerate-)ground state 
l 3 0/TU)) Wltn energy £30 
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o(|3i 
-2a 1 3' 



TUV 



13' 



2T(i)/ 



b\\3' 



13; 



3T(1). 



5T(i)/ 



1 



(A2) 



For brevity the site index T has been suppressed and 



a 
b 



/[12t 2 pd + (e-e p + t ± f 



ale 



e« + 



t±)/(2t pd ) 



(A3) 



For e p ps — 3ey, ps — leV, and £j_ ps 0.25ey one gets 
a ps -0.16, 6 re 0.65, and £ 30 = -7eV. 

In the 2-particle space with no double-occupancy of the 
d-levels there are thirteen states, the creation operators 
of which we label as follows 



9' — o' 

z i; — u i;t u 2;j. z 2i 

2' — d f d f 2' 



rft 2' 
u iU u 2/T z 3i 



U 1U U 2U 

2' — (ft n< 



2' 



2' — d^ 2' — d^ r)l 9' — //t „t 
^13/ = Pjt^U 



(A4) 



To simplify matters we will consider the high-energy 
states with two holes on the oxygen site as decoupled 
from the remaining Hilbcrt space. In the following these 
states will be discarded when evaluating the exchange in- 
tegrals. With this simplification the eigenstates |2y) are 
created by the following set of operators 



2i 

2 2 
2 3 
2 4 



2'! 

2 2 
2' 

2' 

1 

1 

71 



(2' 7 -2' 1 



12) 



& - 2' u ) 



2 7 



(25 + 2 6 2 g 



(25 



2' 4- 2' 



2io) 
2io) 



1 

1 

1 

2 

(2' 5 



(2 7 + 2 12 



(2^ 



2' 



(2' 5 + 2' 6 + 2' 9 + 2' 10 ) 
2r — 2q -+- 2i r 



2n 

z 

212 = (2 S - 2g - 2g + 2 10 + 2 2 13 )03 

213 = (2^ - 2' 5 - 2' 10 + 2' 9 + 7 22i 3 )73 



(A5) 



where, as before, the site index has been suppressed and 



ft (71) = ±e P Tti + 
02(72) = 8t pd //3i( 7 i) 

03(73) =/3i (71)/ 



16^ + (e P -tj_) 2 



(A6) 




with upper(lower) signs on the r.h.s of (A6) referring to 
0(7). The eigenenergies are given by 



(A7) 



and have been labeled into ascending order of their nu- 
merical values as relevant to NaV2C>5. 

The 3-particle space is fourfold degenerate with respect 
to H and the eigenstates are created by 



E 2 12 




= E20 


= 2e p 


-0i/2 p 


a -6.95ey 




■ - E 2s 


= E 2 i 


= e p - 


tx ? 


a -3.25eV 


E2g = ■ 


■ = E 2ll 


= i?22 


= e P 4 


tx F 


a -2.75eV 


E 2 13 




= E23 


= 2e p 


+ 7l/2 * 


a -1.8(W 


E 2l =. 


■ = E 2i 


= E24 


= 







4-1 7 — a 1 -! ji-d n j^Pi*i)i i 4o 



(A8) 



^3/ — dllt^llPtfPll ~ ^111^21 iPlfPl l 

where E40 = 2e p . 

To second order in Hi, the exchange integral J is 
obtained from the energy-dependent transverse spin-flip 
matrix-elements J *~(z) and (z) of the correspond- 
ing second-order effective Hamiltonian 



\j- + {z)= (3 ou 3o/ + i T |VF^— W\3 m 3 0l+n ) 

Z Z — Hq 



(A9) 



where W stands for Hi and the energy variable z is zero 
in the evaluation of the exchange integral. The factor 1/2 
in front of J h corresponds to the fact that in H^ Hels > 

Sf^l+i + &i &i+i = 2(^1 ^1+1 + &l ^l+l)- 

To second order the Ra man scattering amplitude is 

obtained by considering Eq. ( A.9) again, however with W 

denoting the current operator, i.e. %je, in this case and 

with z depending on the energy of the incoming/outgoing 

photon Ui n /uJout- Then, 

R~ + ((jJi n , LOout) = J~ + (uj in ) + J~ + (-LU out ) (A10) 

The first term on the r.h.s of the previous equation de- 
scribes the process where first the incoming photon is 
absorbed in going into the intermediate state, while the 
second term describes the process where the intermediate 
state is reach ed by first emitting the outgoing photon. 
Equation (|A9 l) is evaluated using first the transition 



amplitudes Jfi 1 3 i t (4.) 3 0i+i-L(T) ) 



from the 3<X>3-particlc 



ground states into the bare intermediate 2cg)4-particle 
states \fi) as constructed from (A4) and (|A8j) and second 



by projecti ng th e latt er o nto the 2cg)4-particle eigenstates 
of H , i.e. © and © 

J-+(z) = 2 [(3ou3 oi +i T |W|/x)(M|2«4 jW ) 
(2«4 j n. 1 |i/)(i/|W|3oi T 3oi +u 



z - (E 2i + E i} - 2E- 



:>,(]) 



(All) 



5 



The preceding involves 2(4)-particle states on rungs 
1(1 + 1) only, since the hermitian-conjugate exchange path 
involving 4(2)-particle states on rungs 1(1 + 1) can be 
accounted for by the global prefactor of 2, both, for 
W = Hi and W — jb/e. Moreover, since up to a fac - 
tor of (— i)i = 1 the weights in the numerator of ( All ) 
are identical for W — H\ a nd W = jb/e wc consider 
the former only. Table ( |A12| ) lists the bare intermediate 
states and the corresponding weights. 
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12' 4' \ 
I Z 11Z^1Z+1/ 





-a&(*f + 2t pp ) 


18 


12' 4' \ 

l z 12Z^4Z+l/ 


a6(ij + 2t pp ) 





19 


12' 4' \ 
l z i3Z^3Z+i/ 


b 2 tf 
-b 2 t* 




20 


12' 4' \ 
I z i3i 4: 2;+i/ 


b 2 n 



(A12) 



The constraint in ( All ) on the summation over the in- 
dices fi and v reflects the restriction to i nt erm ediate 
s tates with at most one p-hole. Using (A5), (A7), and 
( A12j ) it is a matter of straightforward algebra to show 
that 



A 



-+ 



i=l 



where 



and 



A: 



AEi 


= E 2 o + 


-©40 — 2^30 « 


1.05eV 


AE 2 


= £21 + 


E40 — 2 £30 « 


A.75eV 


AE 3 


= E22 + 


E40 — 2 £30 « 


5.25eV 


AE 4 


= £23 + 


i?40 — 2 £30 « 


6.20eV 


+ _ _ 


6 2 (8W p(i + 3a/3i) 2 ^ 2 






2[16^ 


+ /3i(e p -t x )] 





(A13) 



(A14) 



(A15) 



^2 — 



-0.0246eU 2 - £0.28 WA + 0.803A 2 

cfbh^t] + tpp) 



-0.00504eU 2 - i0.0133el/A - 0.0886A 2 
a 2 b 2 t* 2 

0.000339eV 2 + i0.00388eVA - 0.0111A 2 

6 2 (8top rf -3a 7 i) 2 q 2 
2[m 2 pd - ll (e p -t ± )} 

-0.000181eV 2 - i0.00206eUA + 0.00589A 2 



Note that the ferro/antiferromagnetic signs of the am- 
plitudes at A = are related to the triplet/singlet char- 
acter of the intermediate states. Eg., A^ + corresponds 
to a matrix element where the intermediate states are 
given by 2g, 2io, 2n, all of which are triplets, therefore 
a ferromagnetic sign of A^ + arises. Inserting the nu- 
merical values of A^ into (A13) we get J v (z = 0) « 



0.049ey - 1.492A 2 /eU + i0.542A. 
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